Robustness of the filamentation instability as shock mediator in arbitrarily oriented 

magnetic field 
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The filamentation instability (sometimes also referred to as "Weibel") is a key process in many 
astrophysical scenario. In the Fireball model for Gamma Ray Bursts, this instability is believed to 
mediate coUisionless shock formation from the collision of two plasma shells. It has been known 
for long that a flow aligned magnetic field can completely cancel this instability. We show here 
that in the general case where there is an angle between the field and the flow, the filamentation 
instability can never be stabilized, regardless of the field strength. The presented model analyzes 
the stability of two symmetric counter-streaming cold electron/proton plasma shells. Relativistic 
effects are accounted for, and various exact analytical results are derived. This result guarantees 
the occurrence of the instability in realistic settings fulfilling the cold approximation. 
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I. INTRODUCTION 

Gamma Ray Bursts (GRB's) and High Energy Cosmic 
Rays (HECR's) are two of the most intriguing enigmas 
in contemporary astrophysics [l|. A particulary promis- 
ing scenario explaining both problems consists in the so- 
called "Fireball model" where particles are accelerated 
through Fermi-type shock acceleration in a relativis- 
tic shock generated by collisions between various shells of 
a relativistic ejecta [1, 0]- Highly energetic particles es- 
caping the shock environment ends up as HECR's, while 
those still undergoing Fermi cycles near the shock front 
radiate in the 7 range, explaining the GRB emission 

This scenario has been validated in recent years 
through a series of pioneering Particle-In-Cell (PIC) Sim- 
ulations where the very shock formation, followed by 
Fermi like particle acceleration, has been directly ob- 
served 0, 3- A key part of the proposed mechanism 
is the very internal shock formation from the collision 
of similar density plasma shells, which is believed to be 
mediated by the filamentation instability occurring when 
coUisionless plasma shells are crossing each others 

The filamentation instability studied here, namely un- 
stable modes propagating perpendicularly to the flow 
of two counter-streaming species, is frequently labeled 
"Weibel" in literature 
ist between these two 



, although differences ex- 
When considering a flow- 
aligned magnetic field, this instability has been alterna- 
tively called "filamentation" jLI] or "Weibel" \ld\. Some 
connections could equall y b e drawn with the so-called 
mirror mode instability jl5l , [l6j . The mirror-mode in- 
stability setting has common points with the magnetized 
version of the Weibel instability In both cases, the 
unstable system consists in an anisotropic plasma im- 
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mersed in an external magnetic field. In our case, the 
instability arises from the relative drift of two isotropic 
species in their own rest frame. Weibel and filamenta- 
tion can intertwine when considering relatively drifting 
anisotropic species [l^, [l^ . 

In view of its central role, much attention has been 
devoted in evaluating the robustness of the filamentation 
instability. It has been found for example that tempera- 
ture effects can significantly reduce its growth rate, and 
even cancel it [lOj, illj . In the same perspective, it has 
been determined that a flow aligned magnetic field can 
suppress the instability [20l - l22j . while baryon loading of 
the plasma shells cari_ increase the growth rate even for 
large temperatures [23]. 

The filamentation instability is far from being the only 
unstable mode triggered by counter-streaming plasmas 
[2^ ]. Nevertheless, it is the only kind of mode producing 
electromagnetic turbulence because the rest of the un- 
stable spectrum of such systems is mainly electrostatic 
[25I . [26l |. The occurrence of the instability is thus de- 
sirable if one wishes to witness the growth of magnetic 
fields, even if it does not dominate the unstable spectrum. 

Although the angle 9 between an hypothetical external 
magnetic field Bq and the initial velocity of the colliding 
shells, is likely to be non-zero in a realistic situation, 
numerical or analytical works usually consider 6 = be- 
sides a few exceptions 27, 28]. The goal of this Letter 
is precisely to explore analytically the general case of a 
non flow-aligned magnetic field. In order to simplify the 
analytical treatment, we focus on the role of the filamen- 
tation instability in the internal shock formation phase by 
considering the collision of two identical relativistic cold 
electron/proton plasma shells. A density ratio different 
from unity would be worth exploring, as important con- 
sequences on the growth rate have been reported [20l[29l]. 
Leaving the asymmetric case for later works, the present 
theory can still form a support for many computational 
works where the encounter of identical plasma shells is 
considered j?!. Isl. I30l - l32j ] . The main result, that we think 
is important for the setting considered, is that the mag- 



FIG. 1: (Color online) Growth rate S {up units) in terms of Z 
7 = 10. 



kvi/ujp and Sis = i^t/^p for various 6. The Lorentz factor is 



netic field can cancel the instability only for 9 — 0. For 
6' 7^ 0, the instability growth rate do decline with Bq, but 
reaches a finite limit for Bq — >■ oo. The occurrence of the 
instability is thus assured, at least in the cold limit. 



II. FORMALISM AND DISPERSION 
EQUATION 

We consider the system formed by two infinite iden- 
tical counter-streaming electron/proton plasma shells of 
density n, velocity ±Vfc and Lorentz factor 7 = (1 — 
w^/c^)"^/^. We choose to align the flow with the z axis 
and set the magnetic field Bq in the (x, z) plane with 

(e^jBo) — 9. We neglect in each shell the protons iner- 
tia, given their much higher mass and their Lorentz fac- 
tor equal to that of the electrons. We work from the cold 
fluid and momentum equations for each electron species 
J = 1,2, 
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where = ^jravj is the momentum and electric charge 
of specie j, and q, m the electron charge and mass respec- 
tively. Linearizing these equations and closing the system 
with Maxwell's equations allows to derive the dielectric 
tensor under the form, 



electronic cyclotron frequency. The dispersion equation 
to be solved is obtained equalling to the determinant 
of the tensor T. Figure [T] displays the numerical eval- 
uation of the growth rate in terms of the reduced wave 
vector Z and the magnetic field parameter fls for various 
obliquity angles 9. We recover previous known results for 
the flow-aligned case, as the growth rate saturates to a 
finite value for Z ^ 00, while the entire spectrum is sta- 
bihzcd beyond fls = l3^/2^ At 9 = 7r/4, the large 
Z saturation is maintained, but the magnetic field effect 
is qualitatively different: although the growth-rate drops 
below a critical magnetization, it does not drop to zero 
but seems to saturate to a finite value. Finally, the case 
9 = tt/2 hardly displays any effect of the magnetic field. 
We will now see that it is possible to analytically access 
much of the quantities involved here. 



III. ANALYTICAL EXPRESSION FOR THE 
MAXIMUM GROWTH RATE 

As can be checked on Figure [U for all values oOIb the 
growth rate quickly reaches a constant value for large 
Z, and this constant is the maximum growth rate for 
all Z's. Note that this Z-saturation of the growth rate 
is observed, not demonstrated. Such behavior of the fil- 
amentation instability for cold systems, magnetized or 
not, has already been reported [12„ |20|]. We thus focus 
from now on the value of the growth rate for Z = 00. 
The dispersion function, though lengthy, is eventually a 
polynomial of degree d in Z which can be written under 
the form. 



^xx ^xy T^xz 
T ~ I -^xy '^yy -^j/z 

xz yz zz 



(3) 



where z* is the complex conjugate of z, while the tensor 
elements are reported in Appendix |^ in terms of the 
following dimensionless variables. 
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The right-hand-side of this equation shows that for Z — > 
00, P{Z) = is equivalent to Od = 0, since the a^'s 
don't depend on Z and m — d < —1 for < m < d — 1. 
Surprisingly, the coefficient Ud can be factorized exactly 
as Od = Q1Q2Q3 with, 
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Here, = Aiinq^/m is the electronic plasma frequency '32 — (s 7 + 2/3 ){2x 7 2^1 ^sva 9) 



of a single shell, and lo}, = \q\BQ/mc the non-relativistic Q3 
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FIG. 2: (Color Online) Numerical evaluation of the maximum 
growth rate Sm of the filamentation instability in terms of the 
magnetization parameter Hb for various field orientations and 
7 = 10. The dotted lines indicate the values given by Eq. 
pop . The green dotted line is coincident with the plain one 
for e = 7r/2, and Eq. dTOj gives 5™ = for 61 = 0. 



It is obvious that Qs = only yields stable modes. We 
prove in Appendix [B] that Qi only yields stable modes, 
while Q2 always yields one purely growing mode with 
TZuj = 0, unless 9 — 0. It has been known for long that 
the filamentation instability can be canceled by a flow 
aligned magnetic field. We here find that this result holds 
only for such field orientation. The dispersion equation 
Q2 = can be solved exactly, yielding the Z-asymptotic, 
and maximum, growth rate Sm as. 
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7 (4/3^73 + (cos 261 - 27^ cos^ e)) , 

-4(3^n% sin^ 6. 
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For = 0, the known result is retrieved with 6m = 
/3\/2/7. Setting 9 = tt/2 in the equations above yields 
Sm = /3\/2/7 regardless of Hb, showing (as hinted by 
Fig. [1]) that the maximum growth rate in the case of a 
perpendicular magnetic field no longer depends on the 
field amplitude. Figure [2] shows the variations of Sm in 
terms of fls for various field orientations. At 9 — 0, 
the su ppr ession of the instability is recovered for fls > 
fi\/2rf [331 • The absence of any field effect for = 7r/2 is 
confirmed, while the curves at intermediate orientation 
show the system undergoes a transition around a critical 
il^. For ils <C ^*B: the system behaves similarly to the 
9 = case. But for fis ^ fi^j, a steady growth rate is 
reached, independent of the magnetization. 

The following expansion is readily derived for low mag- 
netization. 



(8) 



showing the critical magnetic parameter VL*^ is. 



cos 9 



(9) 



For large field i^s S> the growth rate approaches the 
limit. 



1 



7 + 7^ cot^ei' 
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which is simply the maximum growth rate for 9 = 
times a correction factor accounting for the geometry. In 
relativistic conditions, the critical field parameter defined 
by Eq. ^ is always larger than 1. 



IV. DISCUSSION AND CONCLUSION 

We thus find that within the limits of the cold regime 
(see below), the occurrence of the instability is ensured 
as long as the field is not strictly flow aligned. Shock for- 
mation mediated by the filamentation instability seems 
therefore quite robust. Besides the cold approximation, 
we also neglected the proton inertia. Accounting for the 
later is unlikely to result in a lesser unstable system, as it 
has so far been found that Baryon inertia rather sustains 
the instability, even in the kinetic regime p3|. 

The main limitation of this work lies therefore in the 
cold approximation. To be valid, the transverse velocity 
spread with respect to the flow Au^ must fulfllls [2^ 



> Awi(5: 



(11) 



where k is the typical fllamentation wavelength. This 
inequality simply states that particles can be consid- 
ered as monokinetic during one growth period, from 
the standpoint of the plane wave propagating normally 
to the flow [1^. Taking k ^ c/ujp and consider- 
ing Maxwell- Jiittner distribution functions for which 
Au_l/c- (2fcBr/7TOc2)i/2 [2^ yields. 



knT <C '^SLrnc^ . 



(12) 



For weak magnetization with VLb <C Sm ^ \/^Jl 
simply gives kBT <^ mc^, implying non-relativistic tem- 
peratures. In the strong field limit Hb ^ ^*b, the condi- 
tion kBT <C mc^/ (1 -I- 7^ cot^ 9) can be more stringent. 

The Maxwell- Jiittner distribution does not allow for 
a separate definition of a transverse and a parallel tem- 
perature [34] • Considering such distribution function like 
the waterbag, we have Av± ~ AP±/m^ [3^, yielding a 
limitation on the transverse momentum spread only. 
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1 + 72 C0t2 61, VLb > ^*B- 

(13) 

Note that for any function of the form /(v) = 
/^(ti_L, Ti)/|| (t;|| , Til ), the transverse spread J v±f(Pv 
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FIG. 3: Qualitative sketch of the polynomial curves Qi{X), 
left, and Q2iX), right. 



above are not fulfilled. 
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only depends on T±. As it appears, a direct application 
of the present results depends of the kind of distribution 
function considered. A fully relativistic kinetic theory 
would then be needed to assess the robustness of the in- 
stability for oblique propagation when the inequalities 



Appendix A: Tensor elements 

We report here the elements of the dielectric tensor 
defined in Eq. ([3]), 
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Appendix B: Proof that only Q2 ~ yields unstable 
modes in Eqs. (6) 

It is appropriate to set X = a;2 in Qi and Q2, and to 
study the resulting second order equation in X. For both 
polynomials, the coefficient of is 27^ > 0. They thus 
both represent parabolas tending to +00 for X — >■ ±00. 

For the polynomial Qi{X), we have Qi(0) = 
-7(473 + (l+72)n| + (72-l)f72jCos26l) < 0, and 
Qi{0) — 417^ sin2 > 0. It is also straightforward to 
show that the extremum Qi{Xo), where Q[{Xo) = 0, is 
always strictly negative. The curve representing Qi is 
thus necessarily qualitatively similar to the one pictured 



on Fig. [3l The equation Qi{X) = has therefore two 
strictly positive real roots, so that Qi{x) = has only 
purely real roots. 

Regarding Q2{X), we have Q2(0) = -4:l3^n%sm^9 < 
0. Depending on the sign of Q'2{0), the curve for this 
polynomial can only adopt one of the forms pictured on 
Fig. 131 At any rate, Q2{X) — has two purely real 
roots of opposite signs. The equation Q2{x) = will 
then have two purely imaginary roots ±z5, which shows 
in addition that the filamentation instability is here again 
purely growing with TZuj = 0. The growth rate 5 is nec- 
essarily strictly positive unless Q2(0) — 0, that is = 0. 
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